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AIRFOIL-CONTOUR MODIFICATIONS BASED ON c-CURVE 
^5ETH0D OF CALCULATING PRESSURE DISTRIBUTION 
By Theodore Theodorsen 



SUrmARY 



A method, baaed directly on the so-called €-curve 
method published originally In 1951 in NACA Report 
No. I4.II , is presented for use in making modifications 
to the shape and pressure distribution of a given air- 
foil. In particular, it may be desirable to remove 
excessive irregularities or local peaks in the dis- 
tribution. In this process it m.ay be required that 
certain parameters of the airfoil be kept unchanged; 
for instance, the angle of zero lift, the ideal lift 
coefficient, or the moment coefficient. From an aca- 
demic viewpoint, an altered distribution cannot be 
^'prescribed" because compliance with the requirement of 
maintaining a Laplacian flow field is involved. A 
prescribed distribution can therefore not be obtained 
by Iteration* The process, however adequate, is 
necessarily one of qualitative modifications. Several 
numerical examples illustrating the use of the method 
are given in the appendix. 



INTRODUCTI ON 



In 1951 'the author introduced the so-called c-curve 
method for calculating the pressure distribution on air- 
foils of arbitrary shape. (See reference 1.) The 
merit of this method depends essentially on the fact 
that the resulting integral relation can be solved by a 
rapidly convergent process. In the present paper the 
problem of effecting changes in a given pressure dis- 
tribution is considered. The method is based directly 
on the important velocity formula (XTI) of reference 1 
rewritten as formula in reference id, a later 

report. Reference to these papers is made repeatedly 
herein with subsequent omission of the details con- 
ceming the use of the €-curve method. 
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Consider a given airfoil. Prom the e-curve of 
reference 1, \|/ and € are obtained; and from formula 
(XIT) of reference 1 the pressure, or the equivalent 
velocity, is found.' The velocity may be written in 
terms of the quantities cp, c, \|;, xi/^, a + c^i, and v 
as in formula (39') of reference 2; 



,1, pin (a + cp) + sin (a + e-p) 



where the reader referred to references 1 and 2 for 
the meaning of the various symbols. 

Now consider a slightly altered pressure distribu- 
tion. This new pre,?sure distribution is obviously 
related to a new airfoil contour. It is pertinent to 
remark here that from a purely mathematical viewpoint 
the new distribution cannot be "prescribed^* unless the 
new airfoil contour also is prescribed. In a potential 
flow without singularities there exist? a unique rela- 
tionship between the contour and the' pressure distri- 
bution in the flow field. A pressure distribution 
cannot therefore be prescribed (mathematically) for 
the simple reason that the associated contour must be 
given in order to prescribe it. Thus the problem of 
specifying rigorously a pressure distribution is 
reduced '*ad absurdum," Prom an academ.ic standpoint the 
so-called inverse problem therefore does not exist as such. 

Certain alterations of a qualitative nature may be 
performed in spite of the fact that a pressure change 
cannot be prescribed. It is the purpose of this paper 
to indicate a method by which qualitative alteration 
may be performed. It will be noted that the present 
method of contour modification will serve" the intended 
purpose of the inverse method. 

NATURE OP ALTERATIONS 

It" is useful to observe that several types of 
independent alteration are possible. By reference to 
the velocity formula ( 1) , ^f or instance , a change in . 

will appear mainly in the multiplying factor e^ and.- 
will thus effect an increase or a decrease in the veloc- 
ities everywhere on the c'oh tour". This change results 
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simply in a series of airfoils of different thickness as 
the major effect. It is interesting to observe that 
neither the angle of zero lift nor the ideal angle of 
attack has been changed in this operation. The quanti- 
ties € and ^1^ occurring in the velocity formula are 
considered available for the original airfoil contour by 
the c-curve method of reference 1. 

The effect of a change in the angle of attack a 
is well known and need not be discussed here. In fact, 
the main interest lies in improving the pressure dis- 
tribution at and near the optimum, or ideal, angle of 
attack. In the following discussion, therefore, the 
proposed changes are perform.ed at the ideal angle of 
attack only. In other words, the pressure distribution 
is exam.ined at the ideal angle of attack, tentative 
changes are proposed, and results are compared at the 
ideal angle of attack. The restriction that the angle 
of zero lift rem.ain unchanged may or may not be imposed. 
For airfoil contours of zero moment coefficient, as used 
in helicopter blades, the restriction may be imposed that 
the m.oment coefficient remain zero. In the following 
section the nature of such changes is examined with 
several types of restriction used to fulfill specified 
requirements. Such changes may be performed in the 
pressure distribution subject to any one restriction or 
to a combination of several simultaneous restrictions. 



METHOD OF CHANGING THE e -CURVE 



The € -curve can readily be obtained as €(cp) by the 
method of reference 1. In most cases it is desirable 
to keep the ideal lift coefficient constant in order to 
obtain improvement at the exact value of the lift. In- 
asmuch as the expression for the Ideal lift coefficient 
contains the factor 

\ (^T - 

this restriction is equivalent to maintaining a fixed 
difference between and € , the values of c at 

the nose and at the tail, respectively. The absolute 
values may or may not be kept the same. If both 

and € are kept constant in the process of change, 
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the Ideal lift coefficient, the ideal angle of attack, 
and the angle of zero lift are retained. This change 
Is purely local and extremelj^ restricted in nature; only 
!r;inor changes will submit to this stringent type of 
restraint. In order to make a larger change, the con- 
dition of constant angle of zero lift may be relaxed but 
the requirement of a constant ideal lift coefficient 
retained. 

An important case of alteration is the ,case in 
which the moment coefficient is kept constant. It is 
shown in reference 2 that the moment depends on the two 
lowest harmonics in the e(cp) -curve. By prescribing an 
alteration A€(cp) containing higher harmonics than the 
second the pressure distribution may be altered without 
affecting the moment coefficient. Here, al so, further 
restrictions may or may not be imposed. In general, 
the more restrictions imposed, the more manipulations 
are required to adjust the € -curve. 



How a tentative pressure change is translated into 
a change in the € -curve will now be indicated. The 
€(9)- and c f cp) -curves are assumed to be available from 
the method of reference 1. 

A pressure variation Ap along the contour may be 
tentatively prescribed. Since this exact change is not 
expected anyway, exact relationships involving Ap need 
not be used. It is seen from the velocity formula that 

1 p 

— v^ or its equivalent, the pressure Pg measured from 
the stagnation pressure, is given very nearly as 



where A is a function of position only. With almllar 
accuracy, therefore. 



TENTATIVE PRESSURE CHANGES 
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and, finally. 




Ps 



v/here the integral is to b.e taken over the range in 
which the tentative pressure change is given. Because 
this pressure change is improperly chosen, the value of 



for the whole range in which the change is given will 
not, in general, become zero nor will the area under the 
A €-curve 



become zero as required by the conditions on € given 
in reference 1, It is of paramount importance at this 
point to repeat that the originally prescribed pressure 
is necessarily unattainable, as is shown by the fact that 
the two foregoing integrals are, in general, different 
from zero. It will be noticed, however, in the following 
discussion that the essential "shape" effect may be re- 
tained. The orocess is simply to make the A^-curve 
conform to the given requirements by a suitable adjust- 
ment involving the least possible change in the general 
form of the Ae-curve. This adjustment is made by 
changing the location of the maximum and minimum points 
on the curve or by e^^tending the curve beyond the 
original range. The area under the A^-curve can also 
be made zero by changing the reference or mean value. 
Two basic ' condi tions must therefore be Imposed on the 
€-curve; namely, that the two foregoing integrals be 
zero. Several examples are treated in the appendix. 

Finally, a pressure change for constant moment coef- 
ficient must be considered. It will be seen from ref- 
erence 2 that the m.oment depends on the two lowest 
harmonics in the €(cp)-curve and the value of The 
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process Is as follows; Prescribe a tentative pressure 
change Ap, find the corresponding , adjust to comply 
with the two basic conditions previously mentioned, and 
c'etermine the following four Integrals 



"Jo 



,2lt 

\ ^ t: I sin <p dcp 




By removing 

A^ sin cp 4- cos cp + A2 sin 2cp + B2 cos 2<p 

from the initial A c-f unction, the resulting L€ is made 
free of the two lowest harmonics. 

In general, €,p will be changed slightly by the 
removal of the two lowest harmonics. By adding a third 
harmonic 

A^ sin 5<P + B:ap cos 5cp 

both the magnitude and the slope of c at the trailing 
edge may be left unchanged. This end is attained by ■ 
choosing the proper values of A^ and B^. Thus, in the 

function 

■* 

^1^ = -A^ sin 9 - B-j^ cos 9 - A2 sin 2cp 

- B^ cos 2q> + Aj sin Jqp + B^ cos 5<p " 
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the constants A^^ and B are determined by making 



dcp 1 

for 

CP = TT + P 

where 

To the second order In p, there results for and 



A3 = i(Ai - 2A2) + 1(8b^ ^ 5B2)P + Y^A^ ■ 5A2)P^-^--- 

Thus, the six constants are known and the desired change 
in the Ac-curve is given. 



CONCLUDING REMARKS 



It has been pointed out that the so-called inverse 
problem does not exist in a strict sense of the term, 
because a possible pressure distribution cannot be pre- 
scribed unless the new airfoil contour is actually given. 
An airfoil corresponding to a given pressure distribution, 
therefore, cannot in general be arrived at by an itera- 
tion process or by any other method. It is shown that 
only certain qualitative modifications may be effected. 
Such alterations fall logically into several independent 
groups. Attention is given to localized variations, in 
which not only the thickness factor but also the ideal 
angle of attack and the angle of zero lift are kept con- 
stant, pf interest, also, are the pressure changes 



8 



NACA ARR No. lL).G05 



performed with the restriction that the^-inoment coefficient 
and the angle of zero lift remain unchanged. This case 
is of in.portance for airfoils used in helicopter bladeSt^ 
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APPENDIX 



EXAMPLES 



Five cases are treated as examples of the airfoil- 
•contour-modi float ion method of calculating pressure dis- 
tribution: 

R.A.F. 15 airfoil (see table I); local changes on 
lower surface; angle of zero lift, ideal angle 
of attack, and ideal lift doefficient kept con- 
stant 

R.A.F. 15 airfoil; local changes on lower surface; 
moment coefficient and angle of zero lift kept 
constant 

Airfoil contour generated from € - 0.1 sin (a -14-5^)* 
= 0.1 (see table IT); ideal lift coefficient 
kept constant; angle of zero lift and ideal angle 
of attack changed 

Airfoil contour same as in case III; angle of zero 
lift kept constant 

Airfoil contour same as in case III; restrictions 
same as in case I as an example of too severe 
restrictions 

Case I is based on the R.A.F. I5 airfoil, for which 
figure 1 shows the shape and the pressure distribution. 
The purpose of the intended alteration is to remove the 
wavy line on the bottom surface. The first step as 
Indicated is to draw a tentative pressure distribution. 
In the curve at the top of figure 2 the corresponding 
tentative change in pressure is plotted against the 
tngle cp . The next step is to draw the adjusted curve 
for v/hich the area 



I. 

II. 

Ill . 

IV. 
V. 



^ d(P = 0 



Thus a closed A€-curve is assured for case I; this curve 
is called the adjusted curve and is shown In the center 
of figure 2. Tt is also necessary to make the area 
under the Ac-curve ' 
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Jht dtp = 0 

This can be done readily without altering Lt^ or Ac^p, 
as shown by the line f.or. case I, for which 

:i€ dqp = 0 

The corresponding A\|/-curve is shown at the bottom of 
figure 2. ' The modified airfoil shape and pressure 
distribution for case I are shown in figure !• Note 
that the pressure distribution actually obtained differs 
from the tentative pressure distribution and that no 
change has occurred in the angle of zero lift or in the 
ideal angle of attack. 

Case TI is also based on the R^A^P. 15 airfoil, 

with the requirement Imposed that the moment coefficient 
and the angle of zero llf t , remain constant. The tenta- 
tive pressure distribution is the same as that used in 

case I. The ^-curve adjusted for zero area and the 

^s 

Ac -curve adjusted for zero area are therefore Identical 
with those of case I. In this case it is necessary to 
comply with the requirement that the first and second 
harmonics be removed and that some third harmonic be 
added to retain the value of t^. As shown in the dis- 
cussion, the function 



L^t = -A^ sin cp - cos ^ - sin 2cp - cos 2cp 

• + Aj sin 5cp + Bj cos Jtp 

is to be added, where 

1 I 
^0 
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and 

By adding A^c to Ac of case I, the Ac-curve called 

case II and the corresponding A\)r-curve in figure 2 are 
obtained. The resulting modified airfoil contour and 
pressure distribution are shown in figure 3. This case 
is best suited for maintaining zero moment coefficient 
in airfoil sections used in autogyros and helicopters. 

The following three cases, cases m to V, are based on 
the airfoil section generated from c = 0.1 sin (<P - l^^^), 
= 0.1. The original and the tentative pressure 

distributions are shown in figure i;(a). In figure i+(b), 
Ap/pg is shown plotted against the angle (p. The 
tentative pressure curve is adjusted for zero area as 
before. The corresponding Ac-curve is marked "Case 
TIT" in figure 5 and the corresponding A^-curve is 
marked ''Case III" in figure 6. Thus far three choices 
have been treated: 

(1) The exact shape of curve III is retained by 
changing the zero line, or reference line, for Ac, 
which changes both €^ and c but retains the dif- 
ference and therefore ^the ideal lift coefficient. The 
resulting airfoil contour and pressure distribution for 
case III are shown in figure 7. 
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Ac cos <p dcp 



Ac sin 2cp dcp 



A€ cos 2<p dcp 
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(2) In case IV the area under the Ac -curve cor- 
responding to the tentative -curve has been made zero 

by extending the range affected beyond the nose. In 
this case, only the angle of zero lift is kept constant. 
The result as compared with the original is shown in 
figure 8, 

(5) In case V the restriction Is purposely made 
too severe by specifying that no change shall occur 
either in the ideal lift, the angle of zero lift, or 

the angle of ideal lift. 

Case V in figure 5 becomes distorted in attempting 
to fulfill the zero-area requirement and the final results 
shown in figure 9 are correspondingly unsatisfactory. 
The conclusion from this example is that, although cer- 
tain requirements are desirable or required, it is not 
always possible to obtain a good solution within the 
limitations of such requirements. In such a case an- 
other basic type more suitable to the purpose must be 
selected. 
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TABLE I 



ORDI NATES FOR R.A.F. 15 AIRFOIL AND 

MODIFICATIONS, CASES I AND II 

[stations and ordinates in percent 
of wing chordi 



station 



Upner surface 



Original 



Case I 



Case II 



Lower surface 



Original Case I 



Case II 



0 

1.25 

2.5 
5.0 

7.5 

10 

15 

20 

To 

90 

95 
100 



0 

1. 65 
2.50 
^68 
a. 50 
ii.81 

5.55 
5-52 
5.^0 

.09 

5.85 
3.05 
2.01 
1.38 
0 



I: 



0 

1.69 

5.U5 
5.28 
5.00 

4.55 
5.90 

3.12 
2.04 
1.58 

0 



0 

-.70 
-.96 

-1.19 
-1.25 
-1.25 
-1.05 

-.75 

-.28 
'Al 

"II 



■0 

-.s6 
-.50 
-1.03 

-1.13 
-1.25 
-1.27 
-1.25 

-1.55 
-1.17 

-i.o4 
-.71 
-.75 

-.60 

-.1x8 

-.35 

0 



0 

-.63 
-.82 

-.99 
-1.00 

-1.02 
-.91 
-.81 
-.78 

-.60 

::^| 

-.68 

-'.it 

0 
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TABLE II 



ORDI NATES FOR AIRFOIL CONTOUR GENERATED FROM 
€ = 0.1 Sin ((p - ii5°), ^I'o = 0.1, AND 
MODIFICATIONS, CASES III, IV, AND V 



Station 



Upper surface 



Original 



Case III 



Case IV 



Case V 



Lower surface 



Original 



Case III 



case IV 



Case V 



0 

1.25 
2.5 

7.5 

10 

15 

20 

IS 

70 

80 

90 

95 
100 



0 

2.15 

5.05 

6.06 

.20 
.00 
8.78 

8.52 
7.52 
.87 
.19 
2.52 

1.55 
0 



0 

2.80 
I1.06 
ii. 98 

5.76 
6.92 

7.72 
8.59 

3.67 

8.19 
7.50. 

5.90 

2.U1 
l.k2 
0 



0 

2.15 
.05 

.Ll 
5.1+1 
6.22 

7-52 

9.28 
9.U0 
8.78 

18 

2.i|8 

1.59 
0 



0 

2.07 

.Qk 
.21 
5.12 
5.90 
7.00 

.71 

•5^ 
8.S7 
8.06 
7.08 

.67 

.05 
2.18 
1.19 
0 



0 

-1.52 
-2.07 
.2.66 

-2.95 
-5.10 

-5.19 

-5.02 

■m 

-1.27 

-.66 
-.21^ 
-.02 
-.01+ 
-.12 
0 



0 

-1.11 

-1.7U 

-1.97 
-2.15 

-2.ii5 

-2.65 
-2.85 
-2.85 
-2.80 

-2.i|l 
-1.95 
-1.58 
-.82 
-.50 
0 



0 

-1.07 

-1.39 

-1.65 
-1.76 

-1.85 
.1.97 

-2.08 

-2.17 
-2.17 

-1.9k 
-1.61 
-1.26 
-.86 

-.5U 

"39 

0 



-1.53 
-1.76 

-2.31 

-2.71; 

-5.97 
-5.60 

■3.92 
-5.98 
-3.62 
-2.76 
-2.00 
-i.iii 

-.1^ 

.18 

0 
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F/gure 2. - Pressure , -dc-, and A^'Cur\ms, coses I and II. ffiAF. /s airfoil,) 
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Figure 3 . - Shape and pressure distribution , case II . 

(RXF, IS airfoil.) 
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Fig. 4 
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Figure 4.- Pressure disfribuiions and pressure changes, case fJl. 
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F/gure S. - curves for cases III ,Dl, and V. 
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Figure 6. - A^'cur\/es for cases III, IV, and V. 
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Fig. 8 
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Figure 8 - Shape and pressure disfribufion ^ case IV. 
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Fig. 9 
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Figure 9.- Shape, and pressure distribution ^ cose V. 



